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Linearly-sloped or 'ramp' potentials belong to a class of core-softened models which possess a 
liquid- liquid critical point (LLCP) in addition to the usual liquid-gas critical point. Furthermore 
they exhibit thermodynamic anomalies in the density and compressibility, the nature of which may 
be akin to those occurring in water. Previous simulation studies of ramp potentials have focused on 
just one functional form, for which the LLCP is thermodynamically stable. In this work we construct 
a series of ramp potentials, which interpolate between this previously studied form and a ramp-based 
approximation to the Lennard- Jones (LJ) potential. By means of Monte Carlo simulation, we locate 
the LLCP, the first order high density liquid (HDL)-low density liquid (LDL) coexistence line, and 
the line of density maxima for a selection of potentials in the series. We observe that as the LJ limit 
is approached, the LLCP becomes metastable with respect to freezing into a hexagonal close packed 
crystalline solid. The qualitative nature of the phase behaviour in this regime shows a remarkable 
resemblance to that seen in simulation studies of accurate water models. Specifically, the density of 
the liquid phase exceeds that of the solid; the gradient of the metastable LDL-HDL line is negative 
in the pressure (p)-temperature (T) plane; while the line of density maxima in the p-T plane has 
a shape similar to that seen in water and extends well into the stable liquid region of the phase 
diagram. As such, our results lend weight to the 'second critical point' hypothesis as an explanation 
for the anomalous behaviour of water. 



I. INTRODUCTION 

Multiple liquid phases are a common feature of the 
phase diagrams of multi-component mixtures lj. How- 
ever, there is a growing body of experimental and com- 
putational evidence to indicate that they can also occur 
in single component systems. Examples of such have, 
to date, been found in a number of elemental systems in- 
cluding (inter alia] Sulphur B, EL Ell, Phosporous [H|(|0> 
Hydrogen [3, 0, Il3 an d Selenium [ll| . Additionally, ten- 
tative evidence has recently emerged for the existence of 
liquid-liquid transitions in molecular liquids such as n- 
butanol and triphenyl phosphite 12]. 

Arguably, however, the most intriguing example of a 
molecular system exhibiting signs of a liquid-liquid (LL) 
transition, is water. Here the "second critical point" hy- 
pothesis [l3| proposes that the LL transition is wholly 
metastable with respect to freezing and that the associ- 
ated liquid-liquid critical point (LLCP) is responsible for 
the celebrated thermodynamic anomalies in the density 
and compressibility in the stable and metastable liquid 
region near the freezing boundary. Support for this pro- 
posal comes from molecular dynamics simulations of the 
(generally successful) TIP5P interaction potential [l4| . 
These find a metastable LL transition and associated crit- 
ical point, with a line of density maxima which seems to 
emanate from near the LLCP. Moreover, it has been sug- 
gested that at very low temperature the LL boundary 
evolves into a transition between low density and high 
density glassy phases. Whilst transitions between amor- 
phous phases of different densities have been observed 
experimentally, their relationship to the liquid phases is 
still a matter of some debate (see eg. 0,0). To date, 
however, the LL transition has not been observed in real 
water, possibly because the metastable lifetime of these 



phases is too small to be resolvable experimentally. 

Notwithstanding the progress in identifying and char- 
acterizing LL phase transitions and thermodynamic 
anomalies across a variety of disparate systems, it re- 
mains unclear as to whether these phenomena are plu- 
ralistic in physical origin, or can instead be traced to a 
single common underlying mechanism. Furthermore, the 
connection between LL transitions and thermodynamic 
anomalies seems at present rather tenuous: in some sys- 
tems density maxima have been reported without (as yet ) 
evidence of a LL transition, as is the case in Si02 |17j : 
while in other (indeed in most) systems for which LL 
transitions have been reported, they appear unaccompa- 
nied by anomalies (see eg. for reviews). Only in 
water do the two phenomena seem fairly firmly linked. It 
is therefore of interest to enquire as to whether there exist 
simple model potentials that captures the general quali- 
tative behaviour of a LL transition and thermodynamic 
anomalies, and to elucidate their properties in detail. 

Work in this direction has concentrated on the so- 
called core-softened potentials, originally introduced by 
Hemmer and Stell . The functional form of these po- 
tentials is engineered to favor two distinct interparticle 
separations -thereby providing impetus for a transition 
between two liquids of differing densities. Core-softened 
potentials can usefully be subdivided into two classes: 
"shoulder" potentials in which the hard core exhibits a 
region of negative curvature, and "ramp" potentials in 
which the hard core is softened via a linear slope. To 
date, the majority of work on core softened potential 
has concentrated on the shoulder form. Simulation and 
theory HJ HI HI Ei| show that these do indeed 
(given a favorable choices of potential parameters) pos- 
sess a metastable LLCP. However, to date, no firm ev- 
idence of thermodynamic anomalies has been reported. 
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(Initial indications of anomalies in 2D shoulder models 
|25j were subsequently shown [2|J to be an artifact of the 
quasi-continuous nature of the 2d freezing transition in 
the case when the solid has a density less than the liquid.) 

In contrast to their shouldered counterparts, ramp po- 
tentials are known to exhibit both a LLCP and thermo- 
dynamic anomalies. They therefore appear a potentially 
more fruitful route to determining whether the qualita- 
tive features of the LL transition in systems such as wa- 
ter can be described by a very simple model, as well as 
for investigating the more general aspects of the rela- 
tionship between the LL phase behaviour and thermody- 
namic anomalies. As we shall show in the present paper, 
a ramp model can indeed capture (to a remarkable de- 
gree) the qualitative features of the metastable LL transi- 
tion and density maxima seen in simulations of accurate 
water models. 



II. RAMP POTENTIALS 

The phase behaviour of particles interacting via an 
isotropic pair potentials in which the steep repulsive core 
is softened by a linear ramp, was first considered 35 years 
ago by Hemmer and Stell Their calculations for 

a one-dimensional system revealed a range of parame- 
ters for which two phase transitions occurred, and they 
speculated that the same might be true in high dimen- 
sions. More recently, interest in such potentials has been 
rekindled following computer simulation and mean field 
studies of ramp potentials in two and three dimension 
by Jagla [22j . These revealed evidence of HDL and LDL 
phases in addition to the expected liquid and gas phases, 
and the presence of density and compressibility anoma- 
lies. A subsequent detailed MC simulation study of the 
same system by one of us 2£s], accurately mapped a 
portion of the HDL-LDL phase boundary, the liquid-gas 
boundary and the locus of the lines of density and com- 
pressibility maxima. 

The form of ramp potentials we consider in the present 
work is given by: 



which are chosen such as interpolate between the original 
form and a ramp potential approximation to the LJ po- 
tential. The interpolation simultaneously reduces the ra- 
dius of the potential minimum and the maximum range, 
whilst increasing the potential depth. This is done in 
such a way as to maintain approximate constancy of the 
second virial coefficient, thereby ensuring that the po- 
tentials are comparable in a corresponding states sense 

mill. 

We define our family of ramp potentials as follows. The 
limiting value of the potential at the hard core contact 
is held constant at U(r — > r$ ) = 0.69. Choosing to label 
each member of the family of potentials by the radius 
of the minimum ri, the associated well depth D(ri) is 
given by D{r\) = 1.1578 — 0.5578ri, while the value of 
r 2 is tuned such to maintain the second virial coefficient 
at the value B 2 = 1.52. The resulting values of D(ri) 
and r 2 are listed in tab. and a selection of potentials is 
shown in fig. ^ 
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FIG. 1: A selection from the family of ramp potentials studied 
in this work and listed in tab|l] Also shown for comparison is 
the LJ potential whose well depth has been scaled such that 
the second virial coefficient takes the value B2 = 1.52. 
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where U(r) is measured in units of fc^T. Given a con- 
stant hard-core radius ro and contact value U(r — > ) , 
the form of the potential is determined by the position 
of the minimum r% , the maximum range of the potential 
r 2 , and the magnitude of the well depth D = —U{r\). 
In the original work of refs. [26L Wfy , this potential was 
studied for the parameters values r\ = 1.72 r , r 2 = 
3.0 ro, D = 0.1984. In the present work we study the 
properties of a family of such potentials, the members of 



III. COMPUTATIONAL METHODS 

Monte Carlo simulation studies of the phase behaviour 
of a selection of potentials in the family have been per- 
formed within the constant- NpT ensemble. In the results 
reported below, temperature (T) is measured in units of 
the well depth, particle number density (p) is measured 
in units of rj] and the pressure (p) in units of fesT/rg. 
All simulations were performed for systems of N = 300 
particles. 

The principal aim was to locate the liquid-liquid (LL) 
coexistence line and critical point for each potential and 
to map the line of density maxima. The coexistence 
boundaries were obtained using multicanonical Monte 
Carlo techniques and histogram extrapolation in the well 
established manner 133, l3ll 133, 331. To estimate criti- 
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1.61 


0.259756 


2.6438 


1.60 


0.265333 


2.61721 


1.58 


0.276489 


2.56631 


1.5 


0.321111 


2.38847 


1.4 


0.376889 


2.21068 


1.3 


0.432667 


2.06848 



TABLE I: Forms of the ramp potentials studied in this work 
(cf. fig. 0. The well depth is given by D(ri) = —U(ri) = 
1.11578 — 0.568, while the maximum range is tuned such to 
maintain the second virial coefficient at the value B2 ~ 1.52. 

cal parameters, we have employed a crude version of the 
finite-size scaling (FSS) analysis described in ref. [34| . 
The analysis involves scanning the range of pressure p 
and temperature T until the observed probability distri- 
bution of the fluctuating instantaneous particles density 
matches the independently known universal fixed point 
form appropriate to the Ising universality class in the 
FSS limit. Owing to the relatively small temperature 
at which the liquid-liquid critical point generally occurs, 
the acceptance rate for volume updates in the constant - 
NpT ensemble are very low, resulting in extended cor- 
relation time for the density fluctuations. Consequently 
we were able neither to study a wide range of system 
sizes nor obtain data of sufficient statistical quality to 
permit a more sophisticated FSS analysis. Nevertheless 
it transpires that our estimated uncertainties for the crit- 
ical point parameters are sufficient to resolve clearly the 
trends that occur as the form of the potential is altered. 

Lines of density maxima were mapped by measuring 
the density as a function of temperature along isobars of 
the phase diagram. The task of tracking the line of max- 
ima was again aided by histogram extrapolation tech- 
niques. 



IV. RESULTS 

A. The LDL-HDL transition and the density 
anomaly 

It is natural to enquire, first of all, as to the structural 
differences between the LDL and HDL phases. In fig. \2\ 
we show the form of the radial distribution function g(r) 
at a LL coexistence point for the T\ — 1.68 potential. 
The temperature is some 5% below that of the LLCP. 



From the figure, it is evident that in the LDL phase, the 
majority of first neighbors are located at the potential 
minimum, whereas in the HDL phase, there is a much 
larger number of neighbors at the hard core diameter 
and fewer at the potential minimum. 




FIG. 2: Comparison of the forms of g(r) for coexisting points 
on the LL phase boundary. Parameters are r\ = 1.68, T = 
0.0644 = 0.944T C , p = 0.05021. The density of the HDL phase 
is p = 0.484(1), while that of the LDL phase is p = 0.313(1). 

In refs. j2(| [271 it was shown that the ramp potential 
exhibits a maximum in the density as the temperature 
is lowered isobarically though the LDL phase. An ex- 
ample is shown in fig. Pfra) for the case r\ — 1.72 at 
p = 0.0247. For these parameters the density maximum 
occurs at Tmd = 0.095(5). Also included in this figure is 
the form of the radial distribution function g(r) for three 
temperatures on this isobar: one below, one above, and 
one at Tmd- At the hard core contact value, one sees 
that g(r Q ) is greatest for T = Tmd- Thus the anomalous 
density increase with increasing T is apparently due, in 
part at least, to an increase in the number of particles 
choosing the shorter of the two separation distances and 
settling at the hard core, despite the higher energy cost. 

We have traced the locus of the density maxima in the 
p-T plane for several of the potentials studied. These are 
discussed in connection with the phase diagram in the 
following subsection. 



B. Phase behaviour 

Fig.0]shows the location of the LLCP for a selection of 
values of r\ , together with (in some instances) a portion 
of the associated LL coexistence line. One sees that as r\ 
is decreased, the LLCP shifts to lower temperatures and 
higher pressures. On tracking the LL boundary down in 
temperature from the critical point, we observed sponta- 
neous freezing to a hexagonal close packed (hep) struc- 
ture. This solid has a density lower than that of either 
liquid phase [35| 

The freezing point on the LL boundary shifts to higher 
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FIG. 3: (a) The measured number density as a function 
of temperature for p — 0.0247, for the potential having 
ri = 1.72. (b) The measured form of the radial distribution 
function g(r) for the same potential at three temperatures on 
the p = 0.0247 isobar spanning Tmd = 0.095(5). 



temperatures as r\ is decreased. This, coupled with the 
concomitant decrease in the critical temperature, rapidly 
narrows the temperature range over which the LL tran- 
sition is stable as T\ is reduced. For r\ < 1.62, the crit- 
ical point became metastable with respect to the stable 
hep solid phase. For r% = 1.61 and r\ — 1.60, simu- 
lations initiated in the liquid phase was able to sample 
the near critical point fluctuations for a limited time, 
before eventually spontaneous crystallization occurred. 
Freezing occurred very rapidly in the critical region for 
n < 1.60, preventing accurate estimates of critical point 
parameters or indeed the value of n at which the critical 
point becomes completely unstable rather than simply 
metastable. 

It is interesting to note that as n is decreased, the 
initially positive gradient of the LDL-HDL line in the p — 
T plane reduces in magnitude and changes sign close to 
the point at which the LDL-HDL critical point becomes 
metastable. This trend is quantified in fig. [5] Thus the 
gradient of the metastable LDL-HDL line is negative, as 
has also been reported to be the case for water j3(| . We 
find that the change in the sign of the gradient occurs 
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FIG. 4: The near-critical region of the phase diagram for 
each of the ramp potentials studied. Shown in each instance 
is the estimated location of the LLCP. For certain larger val- 
ues of n in the range studied, a segment of the LDL-HDL 
phase boundary has also been estimated. The point of in- 
tersection of the LL phase boundary with the freezing line is 
shown for potentials in which the LLCP is either metastable 
or only moderately stable with respect to freezing. Error bars 
represent the uncertainties in the critical temperature. Un- 
certainties in the critical pressure, as well as in the location of 
the LL line and the freezing points are comparable with the 
symbol sizes. 



because the enthalpy difference between the two phases 
changes sign, rather than the density difference. 




FIG. 5: Estimates of the near-critical gradient of the LDL- 
HDL coexistence boundary in the p — T plane, for the family 
of potentials shown in^] A representative error bar is shown. 

Fig. HJ superimposes the lines of density maxima on the 
phase diagrams of several of the ramp potentials studied. 
The gradient of these lines changes from negative to pos- 
itive values in the p-T plane as pressure is reduced. We 
observe a strong increase in the temperature of the turn- 
ing point as r\ decreases. It is noteworthy that the shape 
of the line of density maxima is similar to that found in 
MD simulations of TIP5P water Furthermore, for 
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the case r± = 1.60, for which the LDL-HDL critical point 
is significantly buried within the stable solid region, the 
density anomaly is nevertheless observable over a wide 
range of the stable liquid phase. 
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FIG. 6: Lines of density maxima for a selection of the po- 
tentials studied, superimposed upon the phase diagrams of 

fig-El 

At high pressure, the lines of anomalies for the various 
potentials becomes rather flat and appear to approach 
the respective LLCP. At low pressure we find that the 
line is truncated by freezing to a face centered cubic (fee) 
solid structure at slightly negative pressures. For most 
of the potentials studied, the anomalous decrease in den- 
sity continues right up to the stable solid region as T is 
lowered isobarically; there is no subsequent density min- 
imum, i.e. a return to "normal" behaviour, as has been 
recently reported in simulations of the ST2 model of wa- 
ter |23. An exception is the the case of r\ ~ 1.72 at 
small negative pressure. Here the density maximum flat- 
tens and a small minimum appears, followed by a clear 
increase in the density as freezing is approached, as is 
shown in fig. 
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FIG. 7: The measured number density as a function of tem- 
perature at p — —0.005, for the potential having r\ — 1.72. 
The system freezes as T = 0.055. 



We have estimated the locus of the liquid-solid coexis- 
tence boundary in the p-T plane for the potential having 
ri = 1.61, for which the LLCP is metastable; the results 
are presented in fig. [SJa) . Fig. OO^b) shows for both a 
high and low pressure, the time evolution of the simu- 
lation density starting from an initial liquid-like config- 
uration for two temperatures either side of the freezing 
point. In the case of the higher pressure (p — 0.1), the 
system freezes to a hep solid of lower density, while for 
lower pressure (p = 0.001), the solid is fee in structure 
having a density greater than that of the liquid. One 
thus expects that the gradient of the freezing boundary 
in the p — T plane is negative at high pressure and posi- 
tive at low pressure. This is indeed confirmed by fig-EIa): 
within the rather limited accuracy of our measurements, 
the gradient of the freezing boundary appears to change 
sign at p ~ 0.02, suggesting that this marks a triple point 
between hep, fee and liquid phases. We have not at- 
tempted to map the hep-fee coexistence line within the 
solid region, although on cooling the fee structure, we 
find it transforms to hep suggesting the gradient of the 
boundary is positive in the p — T plane. 
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FIG. 8: (a) The liquid-solid coexistence boundary and line of 
anomalies in the p-T plane for r± — 1.61. (b) The time evolu- 
tion of the system density close to the freezing transition, as 
described in the text. The figure shows the freezing to a solid 
of higher (lower) density for low (high) pressures respectively. 
Time is measured in units of Monte Carlo sweeps. 
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Whilst the LDL-HDL transition becomes wholly 
metastablc for r\ < 1.62, the lines of density anomalies 
is nevertheless observable in the stable liquid region for 
this value of n and indeed a considerable range of smaller 
ones. However, since no density anomaly occurs for the 
Lennard-Jones potential, it is pertinent to ask how the 
anomaly disappears as we approach this limit. To an- 
swer this question we have studied the case T\ = 1.3 
(cf. fig. which is much closer to the LJ limit than 
the potentials discussed so far. Here we find that freez- 
ing occurs at much higher temperatures than found for 
our studies of the range r\ — 1.72 — 1.62, no anoma- 
lies are seen and there is no indication of a metastable 
liquid-liquid transition. It thus appears that a rapid in- 
crease in the freezing temperature occurs with decreasing 
r*i (as already hinted at in fig. . As a result the stable 
solid region engulfs the temperature range in which the 
anomalies would otherwise occur. This occurs despite the 
fact that the maximum temperature attained by the line 
of anomalies appears to increase slowly as r\ decreases 
(cf. fig. EJ. We estimate that the anomalies are lost for 

n < 1.4. 

V. DISCUSSION AND CONCLUSIONS 

To summarize, previous simulation work on ramp po- 
tentials [2(| [23] has been confined to the situation in 
which the LLCP occupies the stable fluid region. Here 
the LL phase boundary has a positive gradient in the 
pressure-temperature plane of the phase diagram. We 
have shown that by judicious choice of ramp parame- 
ters, one can render the LLCP metastable with respect 
to freezing to a crystalline solid of density lower than that 
of the liquid. A line of density maxima emanates from 
near the metastable LLCP and extends well into the sta- 
ble fluid region. The line of density maxima bends back 
in the p — T plane as pressure is reduced. Furthermore 



(and in contrast to its stable counterpart), the gradient of 
the metastable LL phase boundary is negative. All these 
features are in qualitative agreement with the results of 
simulations of water, and as such, our results lend sub- 
stantial weight to the 'second critical point' hypothesis 
for water. 

It is probably fair to say that there is currently no 
clear picture regarding the factors controlling (i) the ex- 
istence or otherwise of density anomalies in terms of the 
form of the interparticle potential; and (ii) the detailed 
relationship between any such line of anomalies and the 
LL phse boundary. In situations where a line of den- 
sity maxima exists, this is thought to be a sufficient, but 
not a necessary condition for a LLCP to occur, at least 
for supercooled states [40| . However, it remains unclear 
why shoulder potentials exhibit a LL transition, but no 
density anomaly, while ramp potentials exhibit both. As 
regards the locus of the line, thermodynamic considera- 
tions limit the number of ways in which it can terminate 
HH HHi specifically it must either intersect a spinodal 
or transform smoothly into a line of density minima. For 
the family of ramp potentials studied in the present work, 
the line of density anomalies was always found to ap- 
proach the LLCP at their high pressure end. Indeed the 
same appears to be true for a number of other distinct 
models exhibiting LL transitions |4l| , although there are 
yet other models where the intersection appears to occur 
at a point further down the LL boundary [ifl |42]| . We 
have recently obtained preliminary results for the ramp 
potential which may potentiall shed some light on this 
matter. Specifically we find that if the interation range 
is increased to values greater that those studied here, the 
line of density anomalies detaches from the LLCP; its 
intersection with the LL boundary then occurs at sub- 
critical temperatures and in a region of negative pressure. 
We hope to report on this finding in greater detail in a 
future publication. 
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